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Charge neutralisation procedure based on a Markov chain Monte Carlo algorithm 
is applied to a system of generated hadrons. The algorithm changes the charge 
of a randomly picked particle by shifting it within its isomultiplet. For baryons 
changes in both electric charge and baryon number are applied and the algorithm 
leads to charge and baryon number neutralisation. The procedure can thus be 
used to study the effects of the local charge and baryon number conservation. We 
attempt to study these together with the local momentum conservation and their 
effect on azimuthal correlator observable sensitive to local C and CP violation in 
quark-gluon plasma. 



1 Motivation 

Recently, much attention has been paid to the suggestion that local P and CP vio- 
lation can take place in the hot and dense matter produced in heavy-ion collisions. 
In non-central collisions, this leads to charge separation, i.e. an asymmetry in the 
number of positive and negative particles emitted along the angular momentum 
vector of the collision, the so-called Chiral Magnetic Effect [5] . 

The charge separation is described phcnomcnlogically by adding a sine term to 
the Fourier decomposition of the charged-particle azimuthal distribution [3]: 

dN 

= l+2 Vl cos($ a -* flP )+2v 2 cos(2($ a -* i?P )) + ... + 2a ct sin($ Q -V RP ) + ..., 

(1) 

&a, ^rp are the azimuthal angles of the charged particle and the reaction plane, 
respectively (a, (3 standing for the charge of the particles) 

v\, V2 are coefficients linked to the directed and elliptic flows 

a a describes the charge separation relative to the reaction plane. 

Since (a a ) = (ap) = 0, a different observable must be employed to study the 
effect, using the fact that (a a ap) ^ 0. The observable is the azimuthal correlator 
j a f}, defined as follows: 

la p = (cos($ Q + - 2V RP )) (2) 
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= (cos($ Q - V RP ) cos{$p - $ RP )) - (sin($ Q - V RP ) sin($^ - ^ R ,p)) (3) 
= (viaVip) - {a a ap), (4) 

where $ a , are the azimuthal angles of the two charged particles. 

It was, however, observed that the correlator 7 Q/ 3 may be sensitive to local 
conservation of momentum, transverse momentum and charge. These effects have 
been studied separately in various publications |4|5j . We attempt to study the 
momentum, charge and baryon number conservation and their joint effect on 7 Qj g . 



2 Our approach 

The charge conservation is enforced on the system of generated hadrons by a charge 
neutralisation procedure based on a Markov chain Monte Carlo algorithm. The 
charge of a random meson is shifted within its isomultiplet. The process is il- 
lustrated in Figure [1] on the pion isotriplet. The direction of the alteration is 
determined by another random number. 

For baryons, changes in both electric charge and baryon number are applied. 
We show an example of the procedure in Figure [5] on the nucleon isodoublet. Ei- 
ther a nucleon changes into antinucleon or the proton (antiproton) and neutron 
(antineutron) are interchanged. The nature of the alteration is again determined 
by generating a random number. 

The momentum conservation is ensured by generating the particles with the 
Lorentz Invariant Phase Space generator REGGAE pQ. Since the azimuthal cor- 
relator is non-zero only in non-central collisions with elliptic flow being non-zero, 
we simulate the elliptic flow using two samples of hadrons moving away from each 
other in x-direction, as is indicated in Figure [3] 




Figure 1. Charge neutralisation algorithm for Figure 2. Charge and baryon number neutral- 
pions. isation algorithm for nucleons. 
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Figure 3. Toy model simulation of the elliptic flow. 
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Figure 4. Histogram of the simulated integrated elliptic flow V2- 



3 Results 

The variable which we specifically modelled to obtain a system similar to a non- 
central heavy-ion collision is the elliptic flow i>2- We adjusted the parameters of the 
simulation to find values of i>2 corresponding to experimental ones. The histogram 
of integrated «2 is shown in Figure 

The multiplicity dependence of the azimuthal correlator 7 is shown in Figure [SJ 
The correlator was calculated for all pairs of charged particles, like-sign particles 
(positive and negative presented separately) and unlike-sign particles. The three 
types are shown in Figure [5l No significant difference was found between the values 
of the azimuthal correlator of like-sign and unlike-sign particle pairs. This induces 
that the contribution of joint momentum, electrical charge and baryon number 
conservation to the observable 7 is independent of the combined charge of the 
hadron pair. 

The observed signal from the toy model qualitatively agrees with the measured 
correlation of unlike-sign pairs [6]. We miss, however, completely the correlation 
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Figure 5. Multiplicity dependence of the azimuthal correlator 7, triangles and dashed line cor- 
respond to pairs of particles with opposite charge, circles and solid line correspond to negative 
particles, squares and dotted line correspond to positive particles. 



of like-sign pairs. This may be due to non-inclusion of local charge conservation; 
charge was conserved only globally. We shall continue to investigate the topic. 
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